Structure formation in binary mixtures of surfactants: vesicle opening-up to bicelles 

and octopus-like micelles 
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Micelle formation in binary mixtures of surfactants is studied using a coarse-grained molecular 
simulation. When a vesicle composed of lipid and detergent types of molecules is ruptured, a disk- 
shaped micelle, the bicelle, is typically formed. It is found that cup-shaped vesicles and bicelles 
connected with worm-like micelles are also formed depending on the surfactant ratio and critical 
micelle concentration. The obtained octopus shape of micelles agree with those observed in the 
cryo-TEM images reported in [S. Jain and F. S. Bates, Macromol. 37, 1511 (2004).]. Two types of 
connection structures between the worm-like micelles and the bicelles are revealed. 



I. INTRODUCTION 

Amphophilic molecules, such as lipids and detergents, 
self- assemble various structures in aqueous solutions^—. 
In solutions of a single type of simple surfactant, the 
structure of the aggregates is qualitatively understood 
by the relative size of their hydrophilic partsi^. Cone- 
shaped molecules form spherical or worm-like micelles, 
cylindrical molecules form bilayer membranes, and in- 
verted cone-shaped molecules form inverted hexagonal 
structures, or inverted micelles. Mixtures of different 
types of surfactants show more variety of self-assembled 
structures than single-component systems. The edges 
of bilayer membranes can be stabilized by cone-shaped 
surfactants. Hence, discoidal bilayer membranes called 
bicelles (bilayered micelles) or pores on the membranes 
are formed^—. Recently, an undulating cylinder and a 
bilayer disk connected with several worm-like micelles, 
which shape looks like an octopus or jelly-fish, have 
been observed in binary mixtures of diblock copoly- 
mers by cryo-TEMS. When compared with the struc- 
tures of single-component systems, multi-component self- 
assembled structures are not well understood. 

Mixtures of lipids and detergents are used for solubi- 
lization, reconstruction, and crystallization of membrane 
proteins^. In particular, a bicelle is an excellent medium 
to investigate protein functions^^. The proteins main- 
tain their functionality in the bicellesii, and magnetically 
aligned bicelles are used for NMR studies^. The struc- 
tures of several proteins have been determined by the 
bicelle method^^. 

The dynamics of micelle-to-vesicle transitions have 
been studied experimentally by time-resolved scatting 
techniques (light, neutron, and X-ray)^~— . Typically, 
the surfactants aggregate into disk-like micelles, which 
grow and transform into vesicles when their radii ex- 
ceed a critical sizei^. Since fusion between vesicles is 
very slow, the size of the first-formed vesicles is main- 
tained in the time scale of experiments. Thus, the vesicle 
size is determined by kinetics rather than by thermal 
equilibrium. In the opposite transition from vesicles to 
micelles, various pathways of the solubilization of lipo- 
somes are observed depending on the types of lipids and 



detergents^— : Rhythmic shrinkage, bursting, budding, 
peeling, and inside-out inversion. 

Recently, multi-component bilayer membranes have 
been intensively investigated, since they are believed to 
play important roles in living cells. Several patterns of 
phase-separated domains have been observed on giant 
liposomes^—, and the spontaneous curvature produced 
by proteins or glycolipids induce a local high curvature 
such as in a membrane tube^r— . The domain forma- 
tions and morphological changes of the vesicles can be 
observed by optical microscopy. However, the structure 
dynamics of micelle formation are much more difficult 
to investigate; since their size are lOnm scale, electron 
microscopy has to be used for direct observation. 

When molecular simulations are used, the structure 
formation in molecular scales can be revealed in detail. 
In contrast to the fact that many simulations have been 
performed on domain formation in membrane o 31 i 36 ~ — , 
the structure formation of micelles is much less explored. 
In this paper, we study various structures of binary mix- 
tures using a coarse-grained molecular simulation. 

Although computer technology has grown, the typical 
scale for recent simulations of atomistic models is only 
100 ns dynamics of hundreds of lipid molecules. In or- 
der to simulate bilayer membranes on longer and larger 
scales, various types of coarse-grained molecular mod- 
els have been proposed (see review articles^—). Re- 
cently, the bottom-up approaches were taken to construct 
coarse-grained molecular models, where potential param- 
eters are tuned from atomistic simulations^—. On the 
other hand, we choose a top-down approach to construct 
coarse-grained molecular models, where potentials are 
based on the continuum theory. Here, one of the solvent- 
free molecular models^ is employed. In this model, the 
membrane properties such as the bending rigidity and 
spontaneous curvature of the monolayer can be varied 
over wide ranges. 

In Sec. UH the surfactant model and simulation method 
are described. In Sec. IIII1 the membrane properties 
of the phase boundary are investigated using the pla- 
nar membranes with straight edges or pores. The line 
tension of the phase boundary is estimated in Sec. IHI Al 
The phase segregation around the membrane edges and 
the pore shapes on the phase boundary are described 
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in Sees. IfflBl and ILUCl respectively. In Sec. |lVl the 
opening-up of the two-component vesicles are explored. 
Various structures including octopus- like micelles are ob- 
tained. The summary and discussion are given in Sec. fVl 



where f C ut(r) is a C°° cutoff function^, 
/cutM = { e n Xp{A ° (1 + ^)-T)} f<'«*\ (4) 

^ U v — ^cut J 



II. SIMULATION MODEL AND METHOD 
A. Surfactant Model 

A solvent-free molecular models is used to simulate 
two-component surfactant mixtures. We consider N A 
molecules of type A and Nb molecules of type B. The 
type B molecule is a lipid with a cylindrical shape (the 
spontaneous curvature of the monolayer Co — 0). The 
type A molecule is a model of a detergent or lipid 
molecule with a cylindrical or cone shape with a rela- 
tively large head (Co > 0). 

Each (z-th) molecule has a spherical particle with an 
orientation vector Uj , which represents the direction from 
the hydrophobic to the hydrophilic part. There are two 
points of interaction in the molecule: the center of a 
sphere r| and a hydrophilic point r° = + u^cr. The 
molecules interact with each other via the potential, 
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where = r.; - r 3 , nj = |ry|, r y = r^/ry, and fc B T 
is the thermal energy. Here, type A molecules are con- 
sidered for < i < TVa and type B for TVa < i < N — 
N A + N B . 

The molecules have an excluded volume with a 
diameter cr via the repulsive potential, U Iep (r) — 
exp[— 20(r/cr — 1)], with a cutoff at r = 2.4c The second 
term in Eq. (p} represents the attractive interaction be- 
tween the molecules. The two types of molecules can have 
different attractive strengths: e. 



4 t for < i < N A 



and £j = e att f° r < i < N. An attractive multibody 
potential Z/ a tt (pi) is employed to mimic the "hydropho- 
bic" interaction. The potential £/ a tt(Pi) is given by 

C/ att (p i ) = 0.251n[l + exp{-4(p i -p*)}]-C', (2) 

with C = 0.25 ln{l + exp(4p*)}. The local particle den- 
sity pi is approximately the number of particles rf in the 
sphere with radius r at t- 



Pi = ^/cutOlj-). 



(3) 



with n = 6, A = ln(2){(r cut /r att )" - 1}, r att = 1.9a 
(/cut(^att) = 0.5), and the cutoff radius r cut = 2.4er. The 
density p* in U at t(Pi) is the characteristic density, above 
which the pairwise attractive interaction is smoothly cut 
off. 

The third and fourth terms in Eq. ([1]) are discretized 
versions of the tilt and bending potentials of the tilt 
mode l 53 i 54 , respectively. This type of bending potential 
was also used in other coarse-grained lipid model o 55 ' 56 . 
A smoothly truncated Gaussian function^ is employed 
as the weight function 
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with n = 4, r ga = 1.5cr, and r cc = 3cr. The spontaneous 
curvatures of the monolayer membranes are varied by 



r< l 3 
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parameters C bd and C^ d : 
Cm for < t < N A , and C bd = for iV A < i < N. A 
single-component monolayer membrane has spontaneous 
curvature C ~ {fcb cnd /(fc b ond + fctiit)}Cbd/o-— In this 
paper, we use k hcnd = fc tilt , i.e. C ~ C b d/2cr. 

The last term in Eq. ([T]) represents the repulsion be- 
tween the different types of molecules by a monotonic 
decreasing function: Uab('*) — ^i/cut(?") with n = 1, 
A = 1, and r cut = 2.4cr, and A 1 = exp[er/(r cut - a)] to 
set £/ab(c) = 1. 



B. Simulation Method 

The NVT ensemble (constant number of molecules 
N , volume V, and temperature T) is used with peri- 
odic boundary conditions in a box with side lengths L x , 
L y , and L z . A Langevin thermostat is employed to keep 
the temperature. The motion of the center of the mass 
rf = (r| + r?)/2 and the orientation are given by 
underdamped Langevin equations: 
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where to and 7 denote the mass and the moment of 
inertia of the molecule, respectively. The forces are 
given by fj = —dU/drf and fj r = —dU/dui with 
the perpendicular component a 1 ^ = a — (a • Uj)uj and 
a Lagrange multiplier Al to keep u| = 1. Accord- 
ing to the fluctuation-dissipation theorem, the friction 
coefficients £q and £ r and the Gaussian white noises 
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gp(i) and gl(t) obey the following relations: the aver- 



age (CJ*)) = and the variance {g&fijgf^fa)) 



^QlQ2 



t 2 ), where a 1 ,a 2 £ {x,y,z} 



and /3i,/3 2 £ {G,r}. 

The results are displayed with a length unit of a, an 
energy unit of k B T, and a time unit of To = C,gg 2 /k B T. 
The Langevin equations are integrated using the leapfrog 
algorith m 51 ' 59 with m = Cg t o> I = Cr r 0) Cr = Cg " 2 ! an d 
At = 0.005t . 



Membrane Properties of Single-Component 
Membranes 



The potential-parameter dependence of the mem- 
brane properties of single-component membranes (bend- 
ing rigidity k, area compression modulus, lateral diffusion 
coefficient, flip-flop rate, line tensions r c d gc of membrane 
edge, and r br of branching junction) are investigated in 
detail in our previous paper s 51 ' 58 . In this model, these 
properties can be varied over wide ranges. 

The membrane has a wide range of fluid phases, and 
the fluid-gel transition point can be controlled by p* . The 
area compression modulus K A can be varied by fctiit- The 
flip-flop rate can be varied by fctiit and fcbend- The bend- 
ing rigidity k is linearly dependent on fc t nt and fcbend as 
n/k B T = 2.2(fc ti it + fcbend) + b E {e) at C bd = 0. The line 
tension r e d ge of the membrane edge can be controlled by 
varying e and Cbd- As Cbd increases, the line tension 
Tcdgc decreases. Consequently, the bilayer membranes 
become unstable at r e d gc — 0, and instead worm-like 
micelles are formed. 

In this paper, we use fcbend = fctiit = 4 or 8, p* = 14, 
and Cfr d = 0, so that type B molecules form a fluid mem- 
brane with n/k B T — 26 or 44 for e^ tt = 2, respectively. 
We use e^ tt = 2 except for Sec. MI Al where s® tt is var- 
ied with e£j. t = e^ tt . The monolayer (leaflet) of type B 
membrane has zero spontaneous curvature. For type A 



1.5 or 2, and various C£ d . For 



molecules, we use 

e^ tt = 2, the molecules have a very low critical micelle 
concentration (CMC); No or a few isolated molecules are 
seen in the simulation box. For e^ tt = 1.5, the molecules 
have a low but finite value of CMC ~ 10 -4 <7 -3 , which 
corresponds to O.lmM for a ~ 2nm. Vesicles and micelles 
can change their sizes by absorbing and dissolving type A 
molecules at e£ tt = 1.5. The simulations with e^ tt = 1.5 
and 2 correspond to typical lipid-detergent mixtures and 
lipid-lipid mixtures, respectively. 

The same number of type A and type B molecules 
are used for simulations of the planar membranes in 
Sec. HH In Sees. ITTTAl and ImTJl N A = N B = 51 2 
are used. To investigate membrane pores (Sec. IIIIB[) . 
iV A = N B = 1, 024, L x = 38a, L y = 36a, e*. = = 2, 
and fcbend = 8 are used. For most simulations of vesi- 
cle opening-up, the total number of molecule is fixed 
as N — 4, 000, and the number ratio is varied. In a 
few cases, = A^b = 1,000 are used. For the open- 
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FIG. 1: Line tension Fab of the phase boundary of membrane 
domains in tensionless membranes, (a) Dependence on C^a- 
The dashed and solid lines represent data for fcbend = 4 and 8, 
respectively. The symbols (□, o) and (o, A) represent data for 
£ab = 1 and 2, respectively, at e^tt — £ att = 2. The symbols 
(x) and (*) represent data for eab = 1 and 2, respectively, 
1.5, Satt = 2, and fcbend = 8. (b) Dependence on 



at £ a tt 



£ab at fcbend = 8, e£ tt = £att = 2, and C^d =0 (o), 0.4 (□), 

' at r A - r B 



0.6 (A), or 0.8 (o). (c) Tab for various £att 
fcbend = 8, Cj^j = 0, and eab = 1 (□) or 2 (o). The error bars 
are smaller than the size of the symbols. 



mg, 



k 



bend 



tilt 



and L T 



Li, 



L, 



lOOcr 



(N/V = 0.004/cr 3 for N = 4,000) are used. The time 
unit is estimated as To ~ 0.1/iS from the lateral diffu- 
sion coefficient ~ 10 _8 cm 2 /s for phospholipids 60 . We 
consider that micelles maintaining their shapes for an in- 
terval of 100, 000r (~ 10ms) are in steady (meta)stable 
states. The error bars of data are estimated from three 
simulation runs. 
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FIG. 2: Surface tension 7 dependence of the line tension Tab 
of the phase boundary at eab = 2, e^ tt — eftt = 2, fcbcnd = 8, 
and Ctd = (o), 0.4 (□), or 0.6 (A). The error bars are 
smaller than the size of the symbols. 



III. PLANAR MEMBRANES 

A. Line Tension of Phase Boundary 

First, we estimate the line tension Tab of the boundary 
between the type A and B membrane domains. Two 
domains are prepared on the planar membrane along the 
xy plane, and two lines of the phase boundary are set 
parallel to the x axis. The line tension Tab is given 

k 38.61,62 

Tab - (P yy - P XX )V/2L X , (9) 

where V = L x L y L z and 2L X is the total boundary length. 
The diagonal components of the pressure tensor are given 
by 

QTJ 

P aa = (Nk B T - y ai—)/V, (10) 
^ don 

where a £ {x,y,z}. We calculated Tab at three length 
ratios L y /L x = 1.8, 2, and 2.2; and we checked that 
their differences were less than the statistical error bars. 
The membrane surface tension 7 m b is given by 7 m b = 
— (P yy )L Zl since P zz ~ for solvent-free simulations with 
low CMC. 

The line tension Tab can be controlled by £ab (see 
Fig. [T]). It increases with increasing £ab and decreases 
with increasing C^ d . The C^ d dependence is qualitatively 
similar to that of the line tension r e( j g e of the membrane 
edges. On the other hand, the £ a tt dependence is differ- 
ent; Tab is almost independent of e at t [see Fig. Q^c)], 
while r o< 3g increases with increasing e^tt^ 

The tension Tab slightly decreases with increasing sur- 
face tension 7 m b as shown in Fig. [5J The surface tension 
7mb is varied by changing the projected membrane area in 
the xy plane. The repulsive interaction between different 
types of molecules is likely to be reduced by the slightly 
larger distances of neighboring molecules at larger 7 m b- 



FIG. 3: Snapshots of a bilayer membrane with a pore at e att = 
e? tt = 2, fc bcnd = 8, and C^j = 0.4. (a) £AB = 1. (b) eab = 2. 
The red (yellow) and light blue (light yellow) hemispheres 
represent the hydrophilic (hydrophobic) parts of type A and 
type B molecules, respectively. 
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FIG. 4: (a) Schematic representation of a pore on the phase 
boundary, (b) Pore transverse diameter L polc and (c) the 
ratio L B /La of edge lengths at e^tt = Eatt = 2 , fc bc nd = 8, 
and Cb d = (o), 0.2 (□), 0.4 (A), or 0.6 (o). The solid lines 
with the symbols represent simulation data. The dashed lines 
are obtained from Eq. (|13|l . The error bars are smaller than 
the size of the symbols. 



B. Pore on Phase Boundary 

As the membrane area increases, increasing surface 
tension induces a pore opening in the membrane. In the 
phase-separated membrane, the pore stays on the phase 
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boundary or in the domain with lower r G d ge (see Fig. [3]). 
Here, we investigate the pore position and shape in the 
phase-separated membranes for various C£ d and £ab- 

The pore transverse diameter L p0 re and edge length 
ratio Lq /La of type A and type B membranes are esti- 
mated as follows. First, the pore area is calculated using 
a particle-insertion metho d 52 ' 61 . In the xy plane, 10, 000 
ghost particles with a hard-sphere diameter a are dis- 
tributed randomly, and then the particle placement is ac- 
cepted if no overlap with the projections of the hard cores 
of the membrane particles occurs. Next, the molecules 
closer than 1.5a to any ghost pore particle are considered 
to belong to the membrane edge. The neighbor distance 
between molecules are assumed to be same for the both 
domains, so the length ratio is given by the number ratio 
of molecules on the edge: Lb/ La = Np OIC /Np OIC . Two 
mutual contact points of the domain boundary and mem- 
brane edges are defined as the middle points between the 
nearest neighbor edge molecules of types A and B [sec 
Fig. Ufa)]. The transverse diameter L pore is given by the 
length between the two contact points along the x axis. 
Isolated ghost particles and molecules are not taken into 
account in the above calculations. 

At Cj^ d = 0, the two domains are identical. Thus, the 
average pore shape is symmetric, and the edge lengths 
(La) and (Lb) of type A and type B membranes are 
equal (see Fig. 2]). As C£ d increases, the center of the 
pore moves to the type A domain. As £ab increases, the 
pore elongates along the phase boundary, and (La) / (Lb) 
approaches unity. 

These pore shapes are understood by the Young equa- 
tion and Laplace pressure. Three line tensions are bal- 
anced at the mutual contact points^^ 



Tab = Tacos^a) 
r A sin(0 A ) = r B sin(0 B ) 



T B cos(6> B ), 



(11) 



where Ta and Tb denote the line tensions of the mem- 
brane edges of the type A and type B domains, respec- 
tively [see Fig. Ufa)]. Therefore, the contact angle 9a is 
written as 



cos(6» A ) 



r 2 -i- r 2 

1 AB + 1 A 



r 2 

1 B 



2FaFab 



(12) 



The membrane edges of the domains have constant cur- 
vature radii ta and tb- The edge line tension is bal- 
anced with the two-dimensional Laplace pressure A7: 
Ta = ?"aA7 an d Tb = 7"bA7. The transverse diameter 



and edge lengths of each domain are given by 



: poro = 2r A sin(6> A ) = 2r B sin(6» B ), 
La = 2r A 6> A , L B = 2r B B - 



(13) 



Thus, the pore lengths are obtained from the line ten- 
sions and the surface tension difference A7. Since the 
pore has zero surface tension, A7 is calculated as A7 = 
7mb = ~(Pyy)L z - These analytical results qualitatively 
explain the simulation results (see Fig. 2]). However, 




FIG. 5: Bilayer membrane strip at e att = 2, fcbend = 8, and 
Cbd — 0-4- (a) Snapshot of partially phase-separated mem- 
branes at £ab = 0.4. (b) Probability distribution in the y 
direction, which is taken along the eigenvector of the maxi- 
mum eigenvalue of the 2D gyration tensor. The solid lines 
represent data for the type A at eab = 0, 0.4, 0.5, and 0.6. 
The dashed line represents data for the type B at £ab = 0. 
The error bars are shown at several data points for eab = 0, 
0.4, and 0.5. For eab = 0.6, the distribution calculated from 
a single run is shown. 



L porc and Lb/ La are underestimated in particular for 
the parameter sets close to the wetting transition points. 
These deviations are likely caused by the negligence of 
the thermal fluctuations in the theory or the finite size 
effects in the simulations. 

In order to form bicelles with membrane edges cov- 
ered by type A molecules, the complete wetting condition 
Tb > Tab + Ta is necessary: £ab < 1 at C^ d = 0.4 and 
^ab < 2 at C^ d = 0.6. At Cj^j = 0.8, a pore is unstable 
and worm-like micelles grow from the edge of the type A 
domain. 



C. Edge Effects on Phase Separation 

Next, we consider a membrane strip, which has two 
open edges parallel to the x axis [see Fig. |5fa)]. Since the 
type A domain has low line tension Ta < Tb at C^ d > 0, 
the type A molecules prefer to remain at the membrane 
edges. The probability of type A molecules is higher 
than that of the type B at the membrane edges even for 
£ab = [compare solid and dashed lines for £ab = 
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FIG. 6: Effective line tension V of the membrane edge at 
Effbt — 2. (a) The squares (□) and circles (o) represent data 
for C* d = 0.4 and 0.8, respectively, at fcbcnd = 8. The di- 
amonds (o) and triangles (A) represent data for C* d = 0.4 
and 0.8, respectively, at fcbcnd = 4- The dashed lines with (x) 
represent T = Fa + Tab- (b) The triangles (A) and circles (o) 
represent data for fcbcnd = 4 and 8, respectively, at C^ d = 0. 
The error bars are smaller than the size of the symbols. 




FIG. 7: Sequential snapshots of bicelle formation from a vesi- 
cle at iV A = N B = 2, 000, £att = 1.5, £ab = 1, and Cb d = 0.8. 
(a) t = 0. (b) t = 1, OOOro. (c) t = 1, 200r . (d) t = 4, 000r . 
(e) t = 5, 000r . (f) t = 12,000r . (g) t = 28,000r . The 
snapshots in (a-f) are seen from the same viewpoint. 



in Fig. [SJb)]. The free-energy reduction of the edges 
is balanced with the mixing entropy of the two species. 
Similar concentrating effects were reported in the cases of 
the two-component vesicles^ and a vesicle with a cylin- 
drical micelle^, where components with higher sponta- 
neous curvature are concentrated at the high curvature 
regions. 

The effective line tension T of the membrane edge is 
estimated as 

r = -(P XX )V/2L X . (14) 

For single-component membranes, T gives the line ten- 
sion of the membrane edges. For the phase-separated 
membranes, T gives the sum of the two line tensions: 
T = Ta + Tab [compare dashed and solid lines in Fig. 
GDJa)]. Surprisingly, T decreases with increasing £ab at 
low £ab, where the two surfactants are in the mixing con- 
dition (see Fig. [5]). This decrease is caused by the cover- 
age of the membrane edges by type A molecules, which 
reduces the free energy of the edges. At the phase sepa- 
ration threshold T is minimum, and it increases with in- 
creasing £ab for completely phase-separated membranes. 
At £ab = 0, r is slightly lower than r c d ge for the mean 
spontaneous curvature Cbd = (Cj£j + C b 3 d )/2, because of 
the high concentration of type A molecules at the edges 
[see Fig. Mb)]. 



IV. VESICLE OPENING-UP 

As the spontaneous curvature of the monolayer is 
increased by changing Cbd (Co — Cbd/2cr), single- 
component vesicles are ruptured, and these vesicles sub- 
sequently transform into branched worm-like micelles^. 
For two-component vesicles, the number ratio of the com- 
ponents and the line tension Tab of the phase boundary 
are also key parameters. First, we describe the vesicle 
opening-up for £^ tt = 1.5, where type A molecules have 
finite CMC. The vesicle rupture occurs by changing C£ d 
or Tab- Next, we describe the case for e^ tt = 2, where 
both types of molecules have negligibly small CMC. We 
call the former and latter mixtures lipid-detergent mix- 
tures and lipid-lipid mixtures, respectively. 



A. Lipid-Detergent Mixtures 

Typical examples of bicelle formation in lipid-dctergent 
mixtures at Na/Nb = 1 are shown in Figs. [7] and [51 The 
vesicle is equalibrated at eab = and subsequently £ab is 
suddenly changed to eab = 1 at t = 0. After the change, 
the molecules form domains, and then a pore opens in a 
type A domain [see Figs. Efa), (b)] . The pore grows along 
the type A domain and the vesicle opens up [see Figs. 
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FIG. 8: Time development of the bicelle formation at Na = 
N B = 2,000, ettt = 1.5, £ab = 1, and = 0.7, 0.8, or 1. 
(a) The number of molecules na, (b) asphericity a sp , and (c) 
aplanality a p i of the largest cluster. The same data as in Fig. 
□ are shown for C^ d = 0.8. 



(a) . (b) 




FIG. 9: Snapshots of a cup-shaped vesicle at Na = 400, Nb = 
3,600, eftt = 1-5, eab = 1, and C^ d = 0.8. All molecules and 
type A molecules are shown in (a) and (b), respectively. 



E{c), (d)]. In Fig. [7jd), a bicelle and a worm-like micelle 
are detached from the largest cluster. The type A domain 
within the large bicelle diffuses to the membrane edge and 
fuses the other type A domain surrounding the type B 
bilayer disk [see Fig. Efe)]. As consequence of the domain 
fusion, two worm-like micelles connected with the bicelle 




FIG. 10: Sequential snapshots of bicelle formation from a 
vesicle at iV A = 3, 000, N B = 1, 000, e£ tt = 1.5, cab = 0, and 
Cbd = 1. (a) t = 100r . (b) t = l,000r o . (c) t = 20,000r . 
(d) t — 100, OOOro. All snapshots are seen from the same 
viewpoint. 

are formed [see Fig. [7[f)]. After a while, the worm- 
like micelles are detached, and a large discoidal bicelle is 
formed [see Fig. E£g)]. A movie for another example of 
the opening dynamics at Cbd = 0.8 is provided in ESI 
(Movie 1). In the movie, closure of a small pore occurs 
inside a bilayer disk during the bicelle formation. 

At larger Cbd, the vesicle opening begins earlier in 
average, since a pore can open in a smaller domain. 
The connections between bicelles and worm-like micelles 
have longer lifetimes at smaller Cbd- The time devel- 
opment of the largest clusters is shown in Fig. [H] The 
i-th and j-th molecules are considered to belong to the 
same cluster at r?j < r att . The number of molecules 
n c \ in the largest cluster decreases in a stepwise man- 
ner as the micelles undergo detachment. The shape 
of the clusters can be characterized by the ratios of 
the three eigenvalues Ai,A2,A3 of the gyration tensors 
a a = (l/n c i)J2i{ a i ~ a o)(P! ~ Pg), where the center 
of the mass ac = (l/ n ci)Si a i an d a iP € x,y,z. The 
asphericity a SF r^ and aplanarity a p ^ are defined as 

(Ai -A 2 ) 2 + (A 2 -A 3 ) 2 + (A 3 -A 1 ) 2 

2(A 1+ A 2 + A 3 ) 2 ' U5) 

_ 9AxA 2 A 3 

" pl ~ (A 1 + A 2 + A 3 )(A 1 A 2 + A 2 A 3 + A 3 A 1 )' 1 ' 

to quantify the degrees of deviation from a spherical 
shape and planar shape, respectively. These two quan- 
tities are convenient measures to distinguish spherical, 
discoidal, and cylindrical shapes: (a sp ,a p i) = (0,1) for 
spheres, (a sp , a p \) — (1, 0) for thin rods, and (a sp , a p i) = 
(0.25,0) for thin disks. As a vesicle opens up, a sp and 




FIG. 11: Sequential snapshots of a racket-like bicelle forma- 
tion at N A = N B = 1,000, ettt = 2, £ab = 2, and C£ d = 0.8. 
(a) t = 8,500r . (b) t = 9,000r . (c) t = 10,000r . (d) 
t = 48, 000r . (e) t = 60, 000r . 



a p \ show a rapid increase and decrease, respectively. Af- 
ter the worm-like micelles are detached, the bicelle has a 
disk shape, (a sp ,a p i) ~ (0.25,0). 

When type A molecules are too few to cover the cir- 
cumference of the bilayer disk, a cup-like vesicle can be 
formed. Figure ® shows a stable cup-like vesicle. The 
chemical potential of type A molecules surrounding the 
rim is balanced with the bulk chemical potential. Sim- 
ilar stable cup-like vesicles are experimentally observed 
for liposomes with talins^i 

When vesicles contain many type A molecules, the vesi- 
cles are ruptured even at eab = 0. At Na/Nb = 3, 
the rupture begins as soon as a sufficiently high sponta- 
neous curvature is added by changing C^ d = to 1 (see 
Fig. [T0|) . First, a branched worm-like micelle is formed 
like in the single-component case. Then, however, type 
B molecules assemble into the branching junctions and 
form bilayer membranes unlike in the single-component 
case [see a bicelle with three arms in the top right of Fig. 
HOT c)]. Finally, bicelles and worm-like micelles without 
branches are formed. As discussed in Sec. IIII C[ inhomo- 
geneous curvature inside of micelles induces phase seg- 
regation. Furthermore, the segregation induces morpho- 
logical changes of micelles. Thus, the phase segregation 
and morphological changes can be coupled in micelle for- 
mation. 
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FIG. 12: Time development of (a) asphericity a sp and (b) 
aplanality a p i at Na ~ Nb ~ 1,000, e^ tt = 2, eab = 2, and 
C£ d = 0.8. The same data as in Fig. [TT]are shown. 



B. Lipid-Lipid Mixtures 

In lipid-lipid mixtures, the number of molecules in 
surfactant aggregates are changed principally only via 
fission or fusion of aggregates, since the both surfac- 
tants have negligibly small CMC. Figures [Til and [T2l show 
the opening dynamics from a phase-separated vesicle at 
N = 2,000 and C^ d = 0.8. First, the opened mem- 
brane forms a bicelle. This bicelle maintains its shape 
for 38, 000to. Then, a worm-like micelle grows from 
a bicelle edge, and a racket-like micelle is formed [see 
Figs. [TlTd). (e)]. The excess type A molecules in the 
bicelle rim move to the worm-like micelle. The bicelle 
and racket-like micelle can be distinguished by a sp and 
ct p i (see Fig. [T2")) . The worm-like micelle is not detached 
unlike in the finite CMC simulation at e^ tt = 1.5. We 
checked the racket shape is maintained for 100, 000to. 
At N A = N B = 2, 000 and C£ d = 0.8, the vesicles are 
not ruptured. Since higher curvature of smaller vesicles 
slightly reduces a free-energy barrier of the pore open- 
ing, the vesicle rupture occurs more frequently in smaller 
vesicles. 

For vesicles with a high spontaneous curvature at 
Cj^j = 1, ruptured membranes transform into bicelles 
connected with several worm-like micelles (see Fig. 1131) . 
Movies 2 and 3 in ESI show the formation dynamics of 
the octopus- like micelles in Figs. [TST a) and (c), respec- 
tively. After the change from £ab = to 2, type A 
molecules assemble to form domains. As the domains 
become sufficiently large, the membranes in the type A 
domains are ruptured and the excess molecules form con- 
nected worm-like micelles. These processes are similar to 
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FIG. 13: Snapshots of various obtained morphologies of sur- 
factant aggregates at Na = Nb = 2, 000, Eatt = 2, £ab = 2, 
and Cbd = 1. 



those in lipid-detergent mixtures shown in Fig. [7] How- 
ever, the later processes are different. Type A molecules 
do not form a circular rim surrounding type B domains. 
Instead, worm-like micelles are directly connected to the 
bilayer disk of the type B domains as shown in Fig. I13f a). 
These worm-like micelles are stable and their lengths are 
fixed in long simulation runs. At the side of the connec- 
tion junctions, the membrane edges of the type B domain 
are locally exposed in vapor (solvent) space. 

When multiple pores are opened at the same time, 
small pores are closed and the worm-like micelles on the 
pores are detached. If all the pores are closed, polyhe- 
dral vesicles are formed [see Fig. [TBTb)]. The asymmet- 
ric molecular distribution of the inner and outer leaflets 
can generate a spontaneous curvature of the bilayer in 
the type A domains. Hence, the type A domains re- 
main at the vertices with higher curvatures. The forma- 
tion of similar polyhedral vesicles has been reported in 
experiment o 26 ! 27 and other simulations^^. A large pore 
often remains and a cup-like vesicle connected with sev- 
eral worm-like micelles is formed as a steady state [see 
Fig. [Tffl c) and Movie 3]. 

When the spontaneous curvature of the octopus-like 
micelles is changed from C£ d — 1 to 0.8, type A molecules 



form the rim of the bilayer disk connected with worm- like 
micelles like the micelle shown in Fig. [TTJ The number of 
connected worm-like micelles decreases gradually because 
of the length fluctuations due to thermal diffusion. Type 
A molecules in removed short micelles move to remain- 
ing long micelles, a few long worm-like micelles still ex- 
ist after long simulation runs. The cup-like vesicles with 
worm-like micelles are opened up into flat bicelles. When 
the value of C^ d is turned again set as C^ d = 1, the con- 
nection structure transforms back to that shown in Fig. 
[TBT a). Thus, these are not kinetically trapped structures 
but local equilibrium structures. For high spontaneous 
curvature at C£ d = 1, the saddle structure at the side of 
the connection becomes unstable so that the cylindrical 
structure directly connects to the type B domain despite 
partial exposure of the type B membrane edges. 



V. SUMMARY AND DISCUSSION 

We have found various morphologies of micelles in 
mixtures of two types of amphiphilic molecules. For 
lipid-detergent mixtures, bicelle and worm-like micelles 
without branches are formed in dilute solutions. The 
detergent- type molecules stabilize the edges of the bilayer 
disks. On the other hand, for lipid-lipid mixtures, it is 
found that the worm-like micelles are stably connected 
with bicelles. 

Two types of connection structures between the bicelle 
and worm-like micelles are found. At a medium spon- 
taneous curvature of the type A monolayer, worm-like 
micelles of type A molecules are connected with the cir- 
cular type A domain surrounding the type B domain. At 
a high spontaneous curvature, worm-like micelles are di- 
rectly connected with the type B domain. Worm-like mi- 
celles are connected with each other though the rim con- 
sisting of the type A domain in the first case (the round- 
shaped connection) while the absent for the direct con- 
nection in the second case. These octopus-like morpholo- 
gies resemble the cryo-TEM images of diblock copoly- 
mer mixtures^.. These diblock copolymers have negligibly 
small CMC, similar to phospholipids. Since the connec- 
tion points are round in the cryo-TEM images, the struc- 
ture is likely the former round-type connection struc- 
tures; the rim and octopus arms are made of one type 
of copolymer with a relatively longer hydrophilic chain. 
Our simulation suggests that the octopus-shaped micelles 
are generally formed in binary mixtures of negligibly- low- 
CMC surfactants. Hence, similar octopus-like micelles 
can be constructed in binary mixtures of cone-shaped 
lipids and cylindrical-shaped lipids. 

During the solubilization of vesicles in lipid-detergent 
mixtures, octopus-like micelles appear as temporal struc- 
tures. Such temporal structures were considered in some 
experimental studies?' 2 ^. Our simulation results support 
the existence of these intermediate structures. 

To control micelle morphology, several key quantities 
require to be tuned. The spontaneous curvature of the 
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monolayer or the packing parameter— of the surfactants 
determines stable local structure of each phase domain. 
The ratio of line tensions of the phase boundary and 
membrane edges determines the wetting condition for the 
membrane edges. The contact angles of a pore on the 
phase boundary are given by the Young equation. At 
finite CMC, the excess surfactant can be dissolved from 
the micelles, and the composition of surfactants can be 
adjusted to their stable shapes. On the other hand, at 
negligibly low CMC, the number of surfactant molecules 
can be conserved in micelles. Hence, discoidal and cylin- 
drical structures can coexist in a single micelle, where the 
excess amount of cone-shaped molecules to surround the 
circumference form connected worm-like micelles. The 
micelle size can be controlled by kinetics as well as the 
volume fractions of the surfactants. The lipid domain 
size at the time of vesicle rupture determines the result- 



ing bicelle size. 

In this paper, we focus on micelle formation from a 
vesicle. The size and shape of micelles depend on ini- 
tial conditions. Self-assembly from molecules randomly 
dissolved in solutions produces smaller sizes of micelles. 
Micelle formation by self-assembly or other initial condi- 
tions will be explored in further studies. 
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